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ABSTRACT 


Since  the  detection  of  a seismic  event  depends  on  the  signal- 
to-nuise  ratio  at  a receiving  site  or  station  one  may  derive  detection  thresh- 
old magnitudes  from  the  ambient  noise  levels.  The  statistical  distributions 
of  either  peak  noise  amplitudes  or  RMS  noise  levels  then  determine  the  pro- 
bability of  detecting  an  event  of  a given  magnitude.  The  conditions  for  the 
validity  of  this  method  of  estimating  detection  capabilities,  and  the  para- 
meters involved,  are  specified.  The  analysis  shows  that  results  obtained 
from  this  method  should  differ  little  from  maximum  likelihood  direct  method 
estimates.  Application  of  this  method  to  stations  of  the  Very  Long  Period 
Experiment  (VLPE)  confirms  the  analytical  results.  As  few  as  30  reliable 
daily  noise  samples  yield  good  50%  detection  threshold  estimates. 
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SECTION  I 


INTRODUCTION 


In  monitoring  earthquakes  and  underground  nuclear  explosions 
one  is  interested  in  the  detection  capability  of  a seismic  station,  array  or 
network.  This  detection  capability  is  usually  expressed  as  the  probability 
of  detecting  a seismic  event  of  a given  magnitude  in  a certain  region.  One 
method  to  determine  the  detection  capability  is  to  compute  the  percentage  of 
events  actually  detected  for  each  of  several  magnitudes.  The  maximum 
likelihood  approach  to  this  so-called  direct  method  fits  a Gaussian  probability 
of  detection  curve  to  the  detection  percentage s (Ringdal,  1974).  The  dii 'Ct 
method  requires  a large  population  of  events  and  consequently,  establishing 
the  detection  capability  of  a station,  array  or  network  takes  considerable 


time. 


Another  method  estimates  the  detection  capability  based  on  the 
ambient  noise  levels  at  a station  (Lacoss,  1969;  Harley,  1971;  Harley  and 
Heiting,  1972).  In  this  method  it  is  assumed  that  an  event  signal  can  be  de- 
tected when  its  maximum  amplitude  exceeds  that  of  the  surrounding  noise  by 
a certain  margin;  e.  g.  , 3 dB.  We  may  then  assign  a magnitude  to  the  maxi- 
mum noise  amplitude  (A)  occurring  in  a certain  time  gate,  using  some  per- 
iod (T)  and  some  epicentral  distance  (A).  For  surface  waves  this  magnitude 
is  (Harley,  1972): 


M-  . = log  A/T  + log  A + 1.  12  + C. 
snoi  b 1 


(I-D 


where  log  A + 1.  12  represents  the  distance  correction,  and  is  a de- 
tection criterion  margin,  e.  g.  , = 0.  15  for  a 3 dB  margin.  The 
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probability  of  detecting  an  event  of  given  magnitude  Mssig=  x then  is  as" 
sumed  to  be  given  by  the  probability  that  the  noise  magnitude  is  less  than  or 

equal  to  x : 

P(det  M,  . = x)  = P(M  <x)  . (1-2) 

' asig  noi 

Usually,  collecting  noise  data  (for  instance,  an  ensemble  of 
daily  noise  samples)  is  a much  faster  process  than  gathering  data  from  ac- 
tual events,  and  if  noise  amplitude  distributions  are  stationary  with  time  the 
detection  capability  can  be  established  in  a relatively  short  time  period.  This 
is  particularly  important  in  evaluating  new  or  proposed  seismic  station  sites. 

In  the  above  formulae  some  of  the  parameters  involved  have 
not  been  studied  in  detail,  and  so  their  effects  in  estimating  detection  capa- 
bilities have  not  been  fully  specified.  Also,  Lambert  observed  from  his 
Very  Long  Period  Experiment  (VLPE)  data  that  the  noise  peak-ampi  .cde- 
over-RMS  ratio  at  the  seismometer  output  was  relatively  constant  for  the  in- 
dividual stations  as  well  as  from  station  to  station  (Lambert  et  al.  , 1973), 
and  he  raised  the  possibility  that  50%  detection  thresholds  may  be  estimated 
from  average  RMS  noise  levels. 

The  analytical  and  empirical  investigation  of  the  above  topics 
are  the  subject  of  this  report.  Analytically,  the  probability  of  detection  is 
developed  into  validity  conditions  for  the  application  of  detectability  estima- 
tion methods,  and  differences  between  the  estimation  from  noise  and  direct 
estimation  are  discussed.  Empirically,  the  method  of  estimating  detection 
capabilities  from  noise  assuming  a Gaussian  probability  of  detection  is  ap- 
plied to  stations  of  the  VLPE  network,  and  the  results  are  compared  with 
maximum  likelihood  direct  estimates  for  those  stations.  The  study  has  been 
confined  to  surface  wave  detection  threshold  estimates  derived  from  the 
vertical  Rayleigh  wave  component,  but  could  be  extended  to  include  horizon- 
tal Rayleigh  wave,  Love  wave,  and/or  bodywave  detection  threshold  estima- 
tion. 
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Section  II  of  this  report  contains  the  analytical  examination  of 
the  subject.  The  empirical  evaluation,  including  a description  of  the  avail- 
able data,  a discussion  of  the  noise  data  editing  process,  the  detection  thresh- 
old estimation  procedure,  and  the  discussion  of  the  results,  are  given  in  Sec- 
tion III.  The  study  is  summarized  in  Section  IV.  Reference  material  is 
found  in  Section  V,  and  further  details  used  in  this  study  are  given  in  the 
appendices. 
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SECTION  II 

ANALYTICAL  DEVELOPMENT 


1. 

L! 


i 


i 


In  this  section  the  topics  mentioned  in  the  introduction  are  ex- 
amined analytically.  First,  some  surface-wave  magnitude  relationships  are 
established.  Next,  the  probability  of  detecting  an  event  of  given  reference 
magnitude  is  developed  into  the  probability  that  the  event  magnitude  exceeds 
a certain  noise  magnitude.  The  elements  of  this  noise  magnitude  are  specified 
and  examined,  and  their  effects  on  detection  capability  estimation  are  discus- 
sed. Via  a brief  discussion  of  detection  capability  expressions,  this  leads  to 
conditions  and  rules  for  the  validity  of  the  application  of  detection  capability 
estimation  methods,  and  to  specification  of  their  differences.  It  is  then 
shown  that  RMS  noise  levels  may  be  used  in  combination  with  the  maximum 
amplitude-over-RMS  ratio  to  replace  maximum  noise  amplitudes  in  estimating 
detectabilities  from  noise.  Also,  some  considerations  in  averaging  logarith- 
mic terms  are  given.  The  last  part  of  this  section  summarizes  the  analytical 
results. 


A.  SURFACE-WAVE  MAGNITUDE  RELATIONSHIPS 


For  a clear  understanding  of  the  various  surface-wave  magnitude 
values  involved  in  detection  capability  estimation  these  magnitudes  are  defined 
below,  and  their  mutual  relationships  are  given. 


The  surface-wave  magnitude  for  a station  or  array  is  given  by 
the  formula  (Harley,  1972): 


M = log  — + log  A +1.12 


(II- 1) 


who  re : 


J 


A is  the  maximum  peak-to-peak  event  signal  amplitude  in 
m/u  ground  motion  at  the  station,  i.  e.  , the  instrument 
corrected,  or  seismometer  input  amplitude; 

T is  the  period  of  the  maximum  amplitude  A in  sec; 

A is  the  epicentral  distance  in  degrees. 

In  many  cases  "the"  surface  wave  magnitude  is  understood  to 
be  the  magnitude  derived  from  the  maximum  amplitude  with  its  period  near 
20  seconds.  Actual  detections,  however,  may  be  established  from  amplitudes 
with  a different  period.  The-  empirically  determined  magnitude  differences 
due  to  period  effects  are  rather  constant,  so  we  will  use  the  relationship: 

M (T)  = M (2  0)  - d(T)  (H-2) 

s s 

where : 

M (T)  is  the  T-see  magnitude 
s 

M (20)  is  the  20-sec  station  magnitude 
s 

d(T)  is  the  station  magnitude  difference  due  to  period  T. 

Marshall  and  Basham  (1972)  found  a 30-sec  magnitude  differ- 
ence of  0.  30  for  Eurasian  continental  paths.  Lambert  and  Becker  (1973)  com- 
puted a least  mean  square  d(30)  value  for  Northern  hemisphere,  Eurasian 
continental  path  V LPE  stations,  from  both  winter,  and  summer  events;  the 
inter  value  was  0.  3 3,  the  summer  value  0.29,  both  in  good  agreement  with 
Marshall  and  Basham. 

Due  to  radiation  pattern  and  propagation  path  effects  each  sta- 
i i]  .ii  h i ty  will  measure  a different  magnitude  for  one  given  event.  A J -sec 
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network  surface-wave  magnitude  is  obtained  by  averaging,  for  a given  T, 
the  T-sec  magnitudes  measured  at  the  stations  and  arrays  of  the  network 
concerned: 

N 

M*nci(T)  = TT  12  M = i(T) 
i= ! 

where  i is  the  station  index  and  N is  the  number  of  network  stations.  A 
network  may  consist  of,  for  instance,  ail  VLPE  stations  (the  VLPE  network), 
or  may  include  the  LASA,  NORSAR,  or  ALPA  arrays. 

To  evaluate  the  performance  of  a station  or  network  one  must 
have  an  independent,  preferably  global,  reference  magnitude  available.  This 
reference  magnitude  could  be  established  in  the  same  manner  as  the  network 
magnitude,  from  world  wide  stations  and  arrays.  However,  no  current  global 
network  reports  sufficient  M data  to  justify  averaging.  A reference  surface 
wave  magnitude,  therefore,  may  be  obtained  by  converting  a global  bodywaw 
magnitude  to  a 20-sec  surface-wave  magnitude.  For  the  V LPE  network,  for 
example,  the  following  conversion  formula  was  established  by  least -mean 
square  fitting  of  reference  bodywave  magnitudes  and  VLPE  station  20-sec. 
surface  wave  magnitudes  (Lambert  and  Becker,  1973): 

^sref(20)  = 1,20  m|}  ‘ ,-74  (1J-4) 

where  m^  represents  the  reference  bodywave*  magnitude. 

The  station  magnitude  deviation  from  the  20-sec  network  or 
reference  magnitude  is  denoted  by  b in  the  expression 

Ms(20)  - Ms(20)  + b 

in  which  M (20)  and  I)  are  random  variables.  The  mean  deviation  b is  called 
s 

the  station  magnitude  bias. 
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Combining  the  relationships  (II-2 ) and  (11-5)  we  obtain: 


M (T)  = M (20)  - d(T)  + b (11-6) 

S y 


where : 


M (T)  is  the  T-sec  station  magnitude; 
s 

M (20)  is  the  20-sec  reference  magnitude; 
s 

d(T)  is  the  magnitude  difference  due  to  period  *; 

h is  the  station  magnitude  deviation. 

the  prop  ability  of  detection  as  a function  of  event 

MAGNITU  OF. 

One  of  the  criteria  used  by  the  analyst  in  the  mechanism  of 
detecting  seismic  events  is  that  an  event  is  detected  when  at  the  seismometer 
output  the  maximum  signal  amplitude  exceeds  the  maximum  amplitude  of 
the  surrounding  noise  by  a certain  margin  (e.  g.  , * dB).  Other  criteria,  such 

as  dispersion  and  arrival  time  may  also  play  a role  in  the  analyst's  decision 
(e.  g.  , Lambert  ct  al.  , 1973).  To  formulate  the  detection  capability  of  a 

seismic  station,  we  focus  on  the  first  mentioned  criterion,  which  we  will 
call  the  signal -to -noise  ratio  (SNR)  criterion.  The  probability  of  detecting 
events  then  is  determined  by  the  probability  distribution  of  maximum  noise 
amplitudes.  The  probability  of  detecting  an  event  given  its  20-sec  reference 
surface-wave  magnitude  may  be  expressed  as  the  probability  that  &iven  the 
event  magnitude,  its  signal  amplitude  exceeds  the  maximum  noise  amplitude: 


P(drt  M (20)  = x)  - P(ASou|.  > C2-  A^) 


(11-7) 


K 


M (20)  is  the  20-scc  reference  magnitude  for  the  event; 
s ' 

k is  the  given  reference  magnitude  value; 


Ii-d 


ASout 


is  the  maximum  peak-to-peak  event  signal  amplitude 
at  the  seismometer  output; 


is  the  maximum  peak-to-peak  noise  amplitude  at  the 
seismometer  output  within  some  timogate; 


C 


2 


is  the  detection  margin  factor,  e.  g.  , = 1.4  for  a 

3 dB  margin. 


However,  we  would  like  to  express  the  probability  of  detection  as  a probability 
of  magnitude  inequalities.  Since  the  logarithmic  function  is  monotonically  in- 
creasing we  may  replace  the  amplitudes  by  their  base  ten  logarithms  and  in- 
troduce the  parameters  signal  period  (T)  and  epicentral  distance  (A)  to 
establish  a relationship  reflecting  the  event's  station  magnitude.  Since  the 
magnitude  must  represent  the  true  ground  motion  at  the  station,  we  first  cor- 
rect for  the  instrument  response: 


where: 


Ac 

bout 

G(T) 


(H-8) 


is  the  peak-to-peak  event  signal  amplitude  at  the  seis- 
mometer input  in  mJX  true  ground  motion  corresponding 
to  the  maximum  peak-to-peak  output  amplitude; 


T is  the  period  of  the  maximum  output  amplitude; 
G(1  ) is  the  instrument  response  for  the  period  T . 


With  the  modifications  mentioned  above  we  then  obtain: 


F(det  M (20)  = x)  = P(log 

s i 


in 


+ log  A + 1.  12 


- log 


AN 


out 


T*G(T) 


+ log  A + 1.  12  + C ) (II- 9) 


where  C = log  C 
1 < 
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The  left-hand  side  in  the  inequality  now  represents  the  T-sec  station  surface- 
wave  magnitude,  M (T).  Similarly,  we  may  denote  the  right-hand  side  by  a 
magnitude,  and  we  will  call  this  the  station  "noise  magnitude",  : 


A 


N 


out 


M..  - log 

N h T-G(T) 


+ log  A + C 


(11-10) 


where  C = 1.12  + C^.  The  probability  of  detection  then  becomes 


P(det  M (20)  = x)  = P(M(T)>M  1.  (11-11) 

s s N 

Note  that  the  noise  magnitude  is  not  a magnitude  in  the  usual 
sense  because  the  parameters  T,  G(T)  and  A refer  to  the  event  rather 
than  to  the  noise,  and  because  the  noise  magnitude  contains  a detection  mar- 
gin term. 

Using  the  relationship  between  the  T-sec  station  magnitude  and 
the  20-sec  reference  magnitude  (II-6 ) we  obtain: 


P(det  M (20)  = x)  = P(!vl  <20)  - d(T)  + b>  Ml.  (11-12) 

s s [\ 

We  may  now  use  the  substitution  M (20)  - x and  re-arrange  the  terms: 

1 s 

P(det  M (20)  = x)  = P(Ma.  d(T)  - b < x).  (11-13) 

s F 

Thus,  the  probability  of  detecting  an  event,  given  its  20-sec 
reference  magnitude  may  be  expressed  as  the  probability  that  the  random 
variable  (r.v.)  M^  f d(T)  - b is  less  than  or  equal  *o  the  given  reference 
magnitude.  The  statistics  of  all  terms  constituting  this  random  variable, 
therefore,  will  determine  the  statistical  model  on  which  detection  capability 
estimation  methods  may  be  based. 

For  convenience  in  he  further  discussion  we  will  use  the 
symbol  M for  the  r.v.  + d(T)  - b . Furthermore,  we  denote  the 
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probability  distribution  function  of  the  r.  v.  M by  and  obtain  the  final 

relationships: 

P(det  M (20)  = x)  = ^(M<x)  = F fx)  (11-14 

S IVL 

where 

M = M + d(T)  - b . (11-15 

C.  THE  RANDOM  VARIABLE  M = M + d(T)  - I) 

We  will  now  investigate  the  statistical  characieristics  of  each 
term  contributing  to  the  r.  v.  M = M^  t-  d(T)  - I)  . For  this  we  use  the  ex- 
pression for  given  in  Equation  (11-10)  and  slightly  re-arrange  its  vari- 

ables: 

M = log  A nj  - log  T - log  G(T)  + lot;  A + C + d(T)  - b.  (II- 

1Nout 

The  statistical  investigation  thus  will  focus  on  the  terms  log  Aj^  , log  1 
log  G(T),  log  A , d(T),  b , and  any  linear  combination  of  these.  The  term 
C is  a constant  and  needs  no  furthei  examination. 

The  investigation  of  the  r.v.  M and  its  constituents  is  given 
in  Appendix  A.  Histograms  of  actual  data  arc'  shown  for  the  parameters 
loo  Am  , log  T,  and  log  A . The  data  for  AM  reflect  noise  samples 
used  in  the  empirical  evaluation  of  this  study;  the  T a id  A data  were  taker 
from  events  establishing  the  maximum  likelihood  direct  detection  capability 
estimates  (Lambert  et  al.  , 1973)  also  used  in  this  (valuation.  Furthermore 
a relationship  between  the  period  dependent  terms  1,  G(T),  and  d(T)  is  de- 
rived, and  the  statistics  of  the  station  magnitude  variation  are  dismissed. 

The  results  show  to  what  extent  each  of  the  terms  involved, 
and  their  linear  combination,  could  lie  considered  a normal  r.v.  . It  is 
concluded  that  if  the  range  of  relative  distances  is  small  (0.7  < A/A  < Lb. 
where  \ is  the  geometric  moan  distance),  and  bused  only  on  visual 


> > 


inspection  of  the  histograms  concerned,  the  r.  v.  M may  be  considered  to 
be  approximately  normally  distributed  for  most  stations.  Thus,  for  a 
narrow  relative  distance  range,  the  probability  distribution  of  detecting 
an  event  of  giver  reference  magnitude  may  be  considered  to  be  approximately 
Gaussian  in  most  i ases. 

Furthermore,  the  terms  log  G(T)  and  d(T)  were  found  to  be 
approximately  linear  with  log  T.  If  the  parameters  T,  A , and  b are  uncor- 
related, then  the  variance  of  the  r.  v.  M is  approximately: 


M 


2 2 2 2 2 
- cr,  + (/3-«-l)  cr  + a + a 

log  A log  T log  A b 


(H-17) 


whe  re 


a 

log  A 
aiog  T 

log  \ 

a 

b 

Of 


P 


is  the  standard  deviation  of  the  r.  v.  M ; 

is  the  standard  deviation  of  log  Ant  ; 

is  the  standard  deviation  of  log  T ; 

is  the  standard  deviation  of  log  A ; 

is  the  standard  deviation  of  the  station  magnitude  bias; 

is  the  slope  of  the  instrument  response  between  20- 
sec  and  30-sec  periods  ; 

is  the  slope  of  the  d(T)  versus  log  T curve  between 
20-sec  and  30-sec  periods. 


Empirical  values  for  the  above  parameters  result  in  a standard  deviation  range 
for  the  r.  v.  M of 


0.  39*(X  <0.  55  (11-18) 

M 

over  the  stations  concerned.  The  actual  value  of  depends  on  the  station's 
ambient  noise,  its  instrument  response,  the  type  of  propagation  path  between 
the  event  epicenters  and  the  station  (continental  or  mixed  continental-oceanic 
paths),  and  the  relative  distance  range.  Evidently,  for  some  stations  the 


‘J 
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value  0 is  higher  than  the  standard  deviation  (CT=  0.4)  arrived  at  by 
Lambert  et  al.  (1973). 


D.  DETECTION  CAPABILITY  EXPRESSIONS 

To  describe  the  detection  capability  of  a seismic  station,  array 
or  network  one  considers  the  reference  magnitude  values  for  which  the  proba- 
bility of  detection  is  a certain  percentage.  These  reference  magnitude  values 
are  called  detection  thresholds.  For  instance,  the  50%  and  90%  detection 
thresholds,  denoted  by  x^  and  xgg»  respectively,  are  given  by 

P(dct  M (20)  = xcJ  = 0.50  (H-19) 

s 50 

and 

P(det  M (20)  = x ) = 0.90.  (11-20) 

s 90 

Equation  (11-14)  translates  this  into: 


P(M-X50)  = FM(X50>  = °-5°  (II-21) 

and 

P(M<x90)  = Fm(x90)  = 0.90  . (I.r-22) 


Thus,  if  the  probability  distribution  function  of  the  r.  v.  M , 

F (x)  , is  known,  either  in  tabulated,  graphical,  or  closed  mathematical 
M 

form,  the  desired  threshold  value  is  found  from  the  argument  of  F^(x).  For 
instance,  if  M is  a Gaussian  r.  v.  , then  (Papoulis,  1965): 


where : 


is  the  mean  of  the  r.  v.  M ; 
M 


O is  the  standard  deviation  of  the  r.  v.  M . 
M 


(11-23) 
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For  this  case,  the  50%  and  90%  detection  thresholds  are: 


XS0 


(11-24) 


and 


”90 


+ 1. 28  (T  . 

M 


(11-25) 


It  should  be  pointed  out,  however,  that  M need  not  be  Gaus- 
sian to  determine  the  50%  detection  threshold  as  the  mean  of  M . It  suffices 
that  the  probability  density  function  of  M be  even  about  its  mean  to  satisfy 
Equation  (11-24). 


E„  DETECTION  CAPABILITY  ESTIMATION 

In  Subsection  11-B  we  described  mathematically  the  SNR  detec- 
tion mechanism  used  by  the  analyst  in  detecting  seismic  events,  and  we  stud- 
ied the  statistics  involved.  We  arrived  at  a general  detection  formula  which 
is  independent  of  the  method  of  detection  capability  estimation  actually  used. 

Using  this  general  detection  formulation,  we  will  now  describe, 
and  compare  analytically,  two  different  methods  of  detection  capability  esti- 
mation : 

• estimation  from  actual  detection  data  (the  so-called  direct 
method); 

• estimation  from  noise. 

To  each  of  these  methods  two  different  approaches  may  be  taken. 

• make  no  statistical  assumptions; 

• assume  a Gaussian  probability  of  detection. 

E.ich  method's  general  discussion  is  followed  by  a description  of  the  two  ap- 

p roaches.  The  conditions  and  restrictions  for  each  method  and  approach  are 
specified  and  the  expected  differences  in  results  from  either  method  or  ap- 
proach are  stated. 


1. 


The  Direct  Method 


r 


a.  General 

The  direct  method  of  detection  capability  estimation  computes 
the  percentage  of  detections  out  of  ail  ensemble  of  events  in  a certain  region 
and  of  a certain  reference  magnitude.  This  percentage  represents  the  prob- 
ability that  an  event  of  a given  reference  magnitude  can  be  detected.  By  ex- 
pressing the  percentage  as  a function  of  reference  magnitude  one  obtains  the 
detectability  curve,  and  the  ^0%,  90%,  etc.,  detection  thresholls  may  be  de- 
termined from  these  data.  The  method  requires  sufficient  statistical  popula- 
tions of  events  of  the  same  reference  magnitude  in  a certain  region,  over  a 
wide  range  of  reference  magnitudes.  Surface-wave  detection  is  usually  based 
on  SNR,  dispersion  and  expected  arrival  time. 

b.  The  direct  method  without  making  statistical  assumptions 

Applied  in  the  manner  described  above,  the  direct  method  does 
not  require  any  statistical  model  or  assumptions,  and,  in  principle,  the  range 
of  relative  epicentral  distances,  \/\q  , is  not  restricted.  A wide  range,  for 
instance  established  by  events  occurring  over  a large  region  at  a relatively 
short  distance  from  the  station,  probably  would  make  the  detec  tab  ilit  y results 
less  meaningful,  however,  in  the  sense  that  nearby  and  far  away  events  with 
the  same  reference  magnitude  appear  to  have  the  same  probability  of  detection. 

On  the  other  hand,  to  meet  the  requirement  of  a sufficiently  large  population 
of  events  one  may  be  forced  to  estimate  the  detection  capability  from  events 
occurring  over  a relatively  large  region. 

c.  The  maximum  likelihood  approach  to  the  direct  method 

An  .alternative  manner  of  establishing  direct  estimates  is  the  maxi- 
mum likelihood  approach  (Ringdal,  1974).  This  approach  fits  a Gaussian  pro- 
bability of  detection  curve  chrough  the  detection  percentage  points  obtained  by 
the  direct  method,  yielding  the  mean  and  standard  deviation  of  some  threshold 
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magnitude  determined  by  the  seismic  noise  level  and  the  detection  algo- 
rithm characteristics.  The  probability  of  detecting  an  event  of  given  refer- 
ence magnitude  is  then  given  by  the  probability  that  the  event  station  magni- 
tude is  greater  than  the  threshold  magnitude,  given  the  event  reference  mag- 
nitude. 


This  is  similar  to  the  expression  (11-13)  where  the  reference 
magnitude  has  been  translated  into  the  station  magnitude  value  and  the  bias 
term  b . The  threshold  magnitude  in  the  maximum  likelihood  method  then 
is  the  same  as  the  noise  magnitude  plus  the  period  dependent  difference  term, 
+ d(T)  , in  Equation  (11-13). 

The  assumption  that  the  probability  of  detection  is  a Gaussian 
distribution  function  is  approxirr  ately  valid  if  the  range  of  the  relative  epi- 
central  distances,  \/Aq  » is  small,  e.  g.  , for 


0.7<A/Aq<1.6  (11-26) 

|log  A/Aq|  < 0.2  (11-27) 


which  may  establish  an  approximately  log-normal  distribution  of  A . For  ex- 
ample, for  a region  with  Aq  = 50°  , A should  be  within  the  range  35°  to  80° 
to  meet  the  above  constraints. 

Although  the  estimation  model  in  this  approach  is  based  on  the 
signal-to-noise  ratio  detection  criterion,  it  also  depends  cn  the  dispersion 
and  arrival  time  criteria  since  the  method  uses  actual  detection  data  deter- 
mined with  the  aid  of  these  criteria. 

2.  Estimating  Detection  Capabilities  From  Noise 

a.  General 

In  Subsection  II-B  we  found,  according  to  Equation  (11-14)  that 
the  probability  of  detecting  an  event  of  given  reference  magnitude  x equals 
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( 


the  probability  that  the  r.  v.  M = MN  + d(T)  - b is  less  than  or  equal  to  x , 

where  W (Equation  II- 10)  is  the  so-called  noise  magnitude  made  up  by 

N 

logarithmic  terms  of  the  maximum  peak-to-peak  noise  amplitude  at  the  seis- 
mometer output,  the  event  signal  period  T , the  instrument  response  G(T) 
for  that  period,  the  event  epicentral  distance  A , and  a certain  SNR  margin. 
The  term  d(T)  is  the  empirical  difference  between  T-sec  and  20-sec  magni- 
tudes, and  b is  the  station  magnitude  variation  about  the  reference  magni- 
tude. 

The  detection  capability  of  a seismic  station,  therefore,  can 
be  determined  immediately  if  all  statistics  of  the  above  parameters  are  known. 
Thus,  we  collect  noise  samples  to  establish  the  distribution  of  the  logarithm 
of  the  maximum  peak-to-peak  noise  amplitude  at  the  seismometer  output. 
From  existing  data  we  determine  the  distribution  of  the  periods  of  the  maxi- 
mum event  signal  amplitudes.  With  the  known  instrument  response  G(T)  , 
and  the  known  20-sec  minus  T-sec  magnitude  difference  d(T)  , we  then  have 
available  the  distributions  of  log  T,  log  G(T),  and  d(T),  respectively.  We 
then  study  the  seismicity  in  the  area  of  interest  to  obtain  the  distribution  of 
the  logarithm  of  epicentral  distances  between  source  and  station,  log  A , and 
may  determine  the  correlation  between  the  parameters  T and  A . The  dis- 
tribution of  the  station  magnitude  variation  about  the  reference  magnitude, 
finally,  also  is  obtained  from  previous  data  and  its  correlation  with  the  para- 
meters T and  A may  be  studied. 

The  above  procedure  can  be  applied  to  existing  stations  as  well 
as  to  new  or  proposed  sites.  In  the  latter  case  of  course,  a seismometer 
must  be  installed  to  measure  the  noise  amplitudes.  Also,  the  new  station's 
signal  period  distribution  may  differ  from  those  of  other  stations,  due  to  the 
source  radiation  pattern  and  differences  in  propagation  paths.  Furthermore, 
the  new  station's  magnitude  variation  with  respect  to  the  reference  magnitude 
cannot  be  established,  but  may  be  statistically  similar  to  those  of  other 
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stations.  The  instrument  response  of  the  new  station  must  be  known  so  that, 
based  on  the  signal  period  distributions  of  other  stations,  an  indication  about 
the  log  Ci(l)  distribution  may  be  obtained. 

As  pointed  out  in  the  introduction  the  method  of  estimating  de- 
tection capabilities  from  noise  is  fast  compared  with  the  direct  method,  since 
we  only  have  to  collect  noise  data  rather  than  event  signal  data.  With  a well 
functioning  instrument,  and  depending  on  the  time  lapse  between  the  seismic- 
ity sensed  by  the  instrument,  one-hour  noise  samples  can  in  general  be  ob- 
tained at  least  once  a day.  The  empirical  evaluation  of  this  method  described 
in  Section  111  will  show  how  many  samples  are  required  to  obtain  an  accurate 
estimate.  It  is  expected  that  a population  of  approximately  30  valid  noise 
samples  would  suffice.  This  would  mean  that  a station's  detection  capability 
could  lie  determined  within  one  to  two  mon'hs  depending  on  the  quality  of  the 
data.  Estimates  may  have  to  be  adjusted  for  seasonal  noise  differences,  and 
for  station  magnitude  bias. 

b.  Estimation  from  noise  without  making  statistical 
assumptions 

Up  to  this  point  we  have  not  yet  made  or  used  any  statistical 
assumptions;  and  in  principle  detection  capabilities  may  be  estimated  by  em- 
pirically establishing  the  probability  distribution  function,  F^(x),  of  the  r.  v. 
M by  combining  the  event  signal  period,  distance,  and  station  bias  statistics 
with  the  (independent)  noise  amplitude  statistics.  The  50%,  90%,  etc.  , de- 
tection thresholds  then  are  found  from  the  relationships  F..(x  ) = 0.50  , 

M 50 

F (x  ) = 0.  90,  etc.  Since,  as  far  as  SNR  detection  is  concerned,  the  above 
M 90 

described  method  follows  exactly  the  detection  mechanism  used  to  obtain  the 
detection  data  for  the  direct  method,  the  detection  capability  estimates  from 
either  method  should  differ  very  little.  The  only  differences  lie  in  the  incor- 
poration of  dispersion  and  arrival  time  in  the  detection  criterion  used  in  ac- 
tual detection,  and  in  possibly  insufficient  populations  of  events  of  the  same 
magnitude  to  determine  the  detection  percentages  in  the  direct  method. 
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c.  Estimation  from  noise  assuming  a Gaussian  detection 
probability 

If  the  distribution  of  the  r.  v.  M is  an  even  function  the  50% 
detection  threshold  is  found  by  averaging  the  r.  v.  M,  or  averaging  its  con- 
stituents and  adding  the  averages.  If  the  r.v.  M is  Gaussian  we  can  also 
find  the  other  thresholds  in  closed  form;  they  are  given  by  the  error  function. 
The  90%  threshold,  for  instance,  is  given  by  Equation  (II-Z5).  Thus,  assum- 
ing a Gaussian  distribution  function  and  determining  ihe  detection  thresholds 
from  the  error  function  is  essentially  little  different  from  the  maximum  like- 
lihood approach  to  the  direct  method  of  estimating  detection  capabilities.  In 
this  case  the  estimate  differences  also  stem  from  incorporating  dispersion 
and  arrival  time  in  the  detection  criterion  and  from  the  fact  that  in  that 
method  the  /Li  and  <7  values  are  obtained  from  a curve  fitting  procedure 
rather  than  from  explicit  calculations. 

F.  THE  USE  OF  RMS  VALUES  IN  ESTIMATING  DF.'l  ECTION  CAPABIL- 
ITIES FROM  NOISE 

The  introduction  of  this  report  mentioned  Lambi  rt’s  suggestion 
that  RMS  values  rather  than  maximum  amplitudes  might  be  used  to  estimate  a 
station's  detection  capability  from  noise.  In  this  subsection  wo  will  investi- 
gate the  validity  of  that  procedure. 

Lambert  et  al.  (197  1)  observed  that  at  tin  seismome  ter  output, 
the  maximum-noise-amplitude-over- RMS  ratio  was  relative!}  constant.  Wo 
may  change  the  expression  for  the  noise  magnitude  (11-10)  into; 

Mn  ' lo^  HMSoll  4 ioK  R\fs 't-  ' l0^  K'(  l)  + N ' C- 

out 

II,  for  the  case  of  in  even  detection  probability  density  function,  these  terms 
arc  averaged  to  obtain  the  50'  detection  threshold,  we  ma\  make  tin  substi 
tution: 
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since  the  ratio  varies  little  (Appendix  C).  Thus,  if  the  average  value  of 

Am  /RMS  is  known  the  50%  detection  threshold  value  is  determined  by 
Nout  out 


lNout  out 

the  RMS  noise  levels: 


*50  = MN  + d(T)  ' b 


(11-30) 


where 


M = log  RMS  + log  , °U-~  - log  T*G(T)  + log  A + C.  (II  31) 

N out  KMo 


This  method  may  be  convenient  if  for  a given  station  RMS  noise  values  rather 
than  maximum  noise  amplitudes  at  the  seismometer  output  are  available,  and 
the  amplitude-over-RMS  ratio  is  known.  Estimating  the  90%,  etc.,  thresh- 
olds is  more  complicated  since  it  requires  calculating  the  variance  and  co- 
variance  of  the  first  two  terms  in  Equation  (11-31). 


G.  ESTIMATING  THE  50%  DETECTION  THRESHOLD  FOR  AN  EVEN 

DETECTION  PROBABILITY  DENSITY  FUNCTION 


In  this  subsection  we  further  consider  the  process  of  estimating 
detection  capabilities  from  noise.  In  particular,  we  focus  on  the  averaging 
process  involved  in  estimating  the  50%  detection  threshold  from  noise  in  the 
case  that  the  detection  probability  density  function  is  even  about  its  mean. 


For  an  even  distribution  of  the  r.v.  M = M + d(T)  - b the  50%  detection 

N 


threshold  is  given  by  the  mean  of  M , Equation  (11-24): 


x = u = Mxt  + d(T)  - b . 
50  ^ M N 


(11-32) 


Substituting  the  expression  for  the  noise  magnitude  M^  , Equation  (11-10),  and 
following  Equation  (11-16)  we  obtain. 
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X50  = l0g  ANout  " log  T " log  + log  ^ 

U + df(T)  - b + C . (11-33) 

( 

Thus,  we  must  first  compute  the  logarithm  of  the  parameters  , T , 

G(T)  , and  \ , then  average  the  logarithms  for  each  of  these  parameters, 
and  algebraically  add  these  averages  to  d(T)  , b , and  C.  Taking  the  aver- 
age of  the  logarithms  is  equivalent  to  taking  the  logarithm  of  the  geometric 
mean.  In  Appendix  C it  is  shown  that  for  a population  with  a small  cr / fl 
ratio,  the  geometric  mean  approximately  equals  the  arithmetic  mean.  The 
geometric  mean  has  the  characteristic  that  it  weighs  extreme  values  less 
heavily. 

l . According  to  the  above,  the  terms  log  1 and  log  G(T)  may 

be  replaced  by  log  T and  log  C (Tq),  respectively,  where  Tq  is  the  geo- 
metric mean  of  the  signal  period  T.  Since,  furthermore,  d(T)  is  approxi- 
mately linear  with  log  T,  the  term  d(T)  may  similarly  be  replaced  by  d(TQ). 
Alternatively,  we  may  express  the  50%  threshold  as  a function  of  T: 

X50(T)  = 108  ANout  ' log  T ' log  G(T)  + log  A 

+ d(T)  - b + C (11-34) 

where  now  the  period  T is  a running  parameter  rather  than  a mean  value. 

Similar  considerations  hold  for  the  epicentral  distance  term. 

According  to  Equation  (11-33)  we  must  take  the  logarithm  of  its  geometric 

mean,  A : 
o 

log  A = log  Aq  . (11-35) 

Since  the  log  A distribution  is  in  general  not  normal  we  want  to  either  con- 
fine the  \/\  range,  or  make  \ a running  parameter: 
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x50(A)  = ‘0|!  Xut 


log  T*G(T)  + log  A 


+ d(T)  - b + C . (11-36) 

Expressing  the  detection  thresholds  as  a function  of  T and  \ creates  fam- 
ilies of  detection  capability  curves,  each  curve  representing  a certain  period 
and  distance. 


H.  SUMMARY  OF  THE  ANALYTICAL  RESULTS 

Based  on  the  SNR  detection  criterion,  the  foregoing  subsections 
developed  the  probability  of  detecting  a seismic  e.ent  of  a given  reference 
magnitude  x into  the  probability  expression: 


whc re : 


P(det  M (20)  = x)  = P(M  < x)  = F (x)  , 

S M 


(11-37) 


M (20)  is  the  event  20 -sec  reference  magnitude; 
s 


M 

M 


N 


= Mn  + d(T)  - b ; 


is  the  "noise  magnitude"  ; 

T is  the  period  of  the  maximum  event  signal  amplitude 

at  the  seismometer  output  ; 

d(T)  is  the  T-sec  magnitude  correction  with  respect  to  the 
20-sec  magnitude  ; 

b is  the  station  magnitude  variation  about  the  reference 

magnitude  ; 

F (x)  is  the  probability  distribution  function  of  M . 


The  "noise  magnitude"  is  given  by: 


MN  ' 106  ANout 


- log  T*G(T)  + log  A + C , 


(11-38) 
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where: 


Am  is  the  maximum  pcak-to-peak  noise  amplitude  at  the 

Nc  ut 

seismometer  output  ; 

G(T)  is  the  instrument  response  for  the  signal  period  T ; 

\ is  the  epicentral  distance  between  the  event  source 

and  the  station  concerned,  in  degrees  ; 

C = 1.  12  + detection  criterion  margin. 


If  the  statistics  of  all  terms  constituting  the  r.  v.  M are  known 
its  probability  distribution  function,  F (x)  . is  also  known,  and  the  50%,  90%, 
etc.  , ‘detection  thresholds,  x^  Q , x , etc.,  may  be  found  from: 


FM%0>  = °-5°  ' 


(11-59) 


"W  ' °-90- 
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If  F (x)  is  an  even  function  about  the  mean  of  M the  50%  detection  thresh- 

M 

old  equals  this  mean: 


x 


50 


(11-41) 


If,  moreover,  the  r.  v.  M is  Gaussian,  the  detection  thresholds  are  found 
from  the  argument  x of  the  Gaussian  probability  distribution  function 


F 


0.  5 + erf 


x -4. 


M 


(II -4  2) 


where  cr  is  the  standard  deviation  of  the  r.  v.  M.  In  particular,  the  50% 

M 

and  90%  detection  thresholds  are  then  given  by 


X50  “ ^ M 


(H-4I) 
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and 


(11-44) 


x = U + 1.28  & xt  • 

xgO  M M 

It  was  shown  from  available  statistics  that  the  r.v.  M 's  ap- 
proximately a linear  combination  of  approximately  normally  distributed  vari- 
ables or  constants,  if  the  relative  epicentral  distance  range  . AM,,..  kept 
small.  Therefore,  also  the  r.v.  M itself  is  approximately  normally  dt.tr.- 
buted  under  the  above  mentioned  distance  constraint.  If.  moreover,  the  para- 
met-rs  T.  A . and  b are  uncorrelated,  the  standard  deviation  of  M can  be 
computed  from  the  standard  deviation  concerning  each  parameter  involved: 


M 


= aUANout+(^“'1,a‘06T+CT‘OeA 


(11-45) 


where 


(}  = rate  of  change  of  d(T)  with  log  T,  20  T 30  , 
a = rate  of  change  of  log  G(T)  with  log  T,  20  _T  -30  . 
I'rom  available  data  for  well  functioning  VEPE  stations  we  found 


0.  39  < °*55 


(11-46) 


The  above  formulation  is  a model  of  the  SNR  detection  mech- 
anism used  by  the  analyst.  This  model  does  not  take  into  account  the  criteria 
dispersion  and  arrival  time.  Based  on  this  formulation  two  methods  of 
estimating  detection  capabilities  and  two  approaches  to  each  method  were 
cussed-  the  direct  method  and  its  maximum  likelihood  approach:  and  the 
method  of  estimating  the  detection  capability  from  noise,  with  and  without  the 
assumption  of  a Gaussian  probability  of  detection.  It  was  found  that  t ere 
should  be  little  difference  between  the  direct  estimates  and  the  estimates  ob- 
tained from  noise.  Estimating  detection  capabilities  from  noise  is  a relattve- 
ly  last  process,  bu,  it  requires  statistical  information  about  the  signal  periods. 
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the  epicentral  distances  and  the  station  magnitude  bias,  which  information 
may  be  available  from  previously  processed  events.  Accurate  estimates 
should  be  obtainable  from  as  few  as  30  valid  noise  samples.  This  method 
may  use  either  maximum  noise  amplitudes  or  RMS  noise  levels;  either  one 
must  be  measured  at  the  seismometer  output.  In  the  latter  case  the 
maximum-amplitude-over-RMS  noise  ratio  at  the  seismometer  output  must 

be  known. 
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SECTION  III 


EMPIRICAL  EVALUATION 


The  foregoing  theory  was  tested  with  data  available  from  pre- 
vious VLPE  processing.  Noise  samples  taken  by  six  VLPE  stations  during 
the  period  November  1972  through  January  1973  were  edited  for  extreme 
values  and  used  in  estimating  the  detectior  capabilities  of  the  stations  con- 
cerned. The  results  were  compared  with  estimates  obtained  by  the  maximum 
likelihood  approach  to  the  direct  method  (Lambert  et  al.  , 1973).  The 
following  subsections  describe  the  available  data,  the  noise  data  processing, 
and  the  procedure  of  estimating  detection  capabilities  from  noise.  The 
empirical  results,  presented  in  tabular  form,  are  discussed  and  compared 
with  the  conclusions  drawn  from  the  analysis  in  the  previous  section. 

The  locations  of  the  stations  of  the  VLPE  network  are  given 
in  Table  III  - 1 and  Figure  III  - 1 . The  stations  TLO,  OGD,  and  FRK  were  not 
in  operation  during  the  above  mentioned  period,  and  the  data  from  EIL  and 
MAT  were  considered  unreliable. 

A.  AVAILABLE  DATA 

In  the  following  we  will  briefly  describe  the  available  noise 
data  and  maximum  likelihood  detection  capability  estimates  used  to  perform 
the  empirical  evaluation. 

1 . Noise  Data 

Ensembles  of  noise  samples  collected  during  the  period  Nov- 
ember 1972  through  January  197  3 at  each  of  (he  VLPE  stations  CTA,  CHG, 
KON,  KIP,  ALQ,  and  ZLP  were  used  to  estimate  detection  capabilities  from 
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TABLE  III- 1 

VERY  LONG  PERIOD  EXPERIMENT  (VLPE) 
STATIONS  AND  LOCATIONS 


Station 


Designator  Latitude  Longitude 


Charters  Towers,  Australia 
Chiang  Mai,  Thailand 
Fairbanks,  Alaska 
Toledo,  Spain 
Eilat,  Israel 
Kongsberg.  Norway 
Ogdensburg.  New  Jersey 
Kipapa,  Hawaii 
Albuquerque,  New  Mexico 
La  Pax,  Bolivia 
Matsushiro,  Japan 


CT  A 
CHG 
FBK 
T IX) 


KON 

OGD 


ALQ 

ZLP 

MAT 


20.  09S 
18.  79N 


39.  86N 
29.  55N 
59.  65N 
41. 07N 


16.  5 OS 


146.  26E 
98.  98E 


64.90N  148.  01  W 


4.  02  W 
34.  95E 
9.  59E 
74. 62 W 


21.42N  158. 02W 

34.  94 N 106.  46  W 


68.  1 3 W 


36.  54 N 138.  12E 


noise.  The  noise  samples  consist  of  digitized  seismograms  taken  over  ap- 
proximately one-hour  periods  when,  according  to  PDE,  NORSAR,  and  LASA 
listings,  no  seismic  activity  was  present.  The  digital  data  expressed  in  com- 
puter count  (c.  c.  ) units  were  bandpassed  from  0.  0 to  0.  14  Hz  via  Fourier 
transforms,  and  normalized  such  that  amplitudes  with  a 40-sec  period 
represent  true  ground  motion  in  millimicrons.  The  relationship  between 
input  (true  ground  motion)  and  output  amplitudes  as  a function  of  period 
then  is  given  by  the  station  instrument  response  (G(T)  with 

G(40)  = 1.0  , (HI- 1) 


and  the  response  at  all  other  periods  is  scaled  relative  to  G(40).  A typical 
instrument  response  curve  is  given  in  Appendix  A (Figure  A-9);  the  instru- 
ment response  curves  for  all  stations  of  the  VLPE  network  are  given  by 
Lambert  et  al.  (197  5).  Table  A- 1 lists  the  20-sec  and  30-sec  instrument 
responses  for  the  period  that  the  data  were  collected.  Further  details  on  col- 
lecting and  processing  VLPE  noise  data  are  given  by  Prahi  (1974). 

2.  Maximum  Likelihood  Detection  Capability  Estimates 

The  detection  capability  estimates  obtained  from  these  noise 
data  was  to  be  compared  with  estimates  found  by  means  of  the  maximum  like- 
lihood approach  to  the  direct  method  applied  to  actual  detection  versus  non- 
detection  of  events  which  occurred  during  the  periods  January  1,  1972  through 
March  20,  1972;  June  1,  1972  through  August  31,  1972;  and  November  1,  1972 
through  December  31  , 1972.  The  maximum  likelihood  detection  threshold  es- 
timation was  expressed  in  terms  of  reported  bodywavc  magnitudes,  and  the 
bodywave  magnitude  30%  detection  thresholds  were  subsequently  converted  to 
20-sec  surface  wave  magnitude  50%  detection  thresholds  using  the  relation- 
ship 

M (2C*  1.20  ni  - 1.74  . (111-2) 

s b 
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The  same  relationship  is  used  to  establish  a reference  surface  wave  magni- 
tude, Equation  (11-4).  This  conversion  formula  v/as  established  by  least 
mean  square  fitting  of  available  m and  data  by  Lambert  and  Becker 
(1973).  Furthermore,  the  maximum  likelihood  90%  detection  thresholds  were 
determined  from 

x = + 1 . Z 8 cr  =x  + 1.28c  (III-3) 

90  90 

in  which  a (T-value  of  0.4  was  used  for  all  stations.  According  to  the  dis- 
cussion in  Subsection  II-C  this  constant  a-value  used  for  all  stations  may 
have  been  a low  estimate  for  some  stations,  since  for  our  Gaussian  detection 
probability  model  we  found  a range 

0.  39  < <7  <0.  55  < HI “4 ) 

M 

depending  on  the  station  noise  amplitude  distribution,  the  propagation  path, 
the  instrument  response,  and  the  range  of  epicentral  distances. 

B.  NOISE  DATA  PROCESSING 

The  noise  data  described  in  the  previous  subsection  needed  to 
be  bandpassed  to  agree  with  the  event  signal  processing  band,  and  edited  to 
remove  samples  taken  during  station  malfunction  or  station  calibration,  and 
samples  containing  spikes  and  unreported  signals.  The  best  way  of  editing 
the  data  is  visual  inspection  of  the  wide-band  noise  traces.  Even  then,  in 
some  cases,  some  ambiguity  remains  as  to  whether  a sample  contains  high 
level  noise  or  one  of  the  anomalies  mentioned  above.  Furthermore,  CAL- 
COMP  plotting  and  the  visual  inspection  of  the  plotted  traces  for  all  samples 
(some  600  traces)  was  considered  too  time  consuming  (especially  since  an 
interactive  display  capability  was  being  developed  for  the  PDP-1S  computer 
at  the  Seismic  Data  Analysis  Center  (SDAC)  and  would  be  available  for  spot- 
checking data  quality).  Therefore,  the  following  automatic,  but  knowingly 
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crude,  method  was  applied  while  also  preparing  the  data  for  detection  thresh- 
old estimation.  To  facilitate  reference  to  the  computer-generated  Tables 
HI _3  through  III -8,  and  E-l  through  E-6,  the  variables  (AMP,  RMS,  AMP/ 
RMS,  etc.)  presented  in  those  tables  are  acco rdingly  given  in  capital  letters 

in  the  following  text. 

The  available  (wide-band)  noise  data  were  subjected  to  a 17- 
44  sec  bandpass  filter.  This  passband  was  also  used  in  standard  processing 
of  VLPE  event  signals.  From  the  bandpassed  noise  data,  the  maximum  peak 
to -peak  amplitude  (AMP)  and  RMS  values  were  determined.  The  Peak-to- 
peak  amplitudes  were  found  by  a search  algorithm  which  also  yields  the  per- 
iod (AMP  PER)  of  the  maximum  noise  amplitude.  The  RMS  value  is  com 

puted  in  the  time  domain  as 


RMS 


(III-5) 


where 

n(i)  is  the  ith  data  poin'  noise  value  in  mji  ; 

N is  the  total  number  of  data  points  in  the  time  trace. 

Next,  the  values  AMP/RMS,  LOG  AMP  , LOG  RMS,  LOG 
RATIO  (=  log  AMP/ RMS)  were  computer  for  each  sample.  The  mean  and 
standard  deviation  of  the  values  AMP,  RMS,  AMP/ RMS,  LOG  AMP,  LOG 


RMS,  and  LOG  RATIO  were  determined  according  to 


and 


MEAN 


N 


d(i) 


STD.DEV. 


( III- 6) 


( HI  - 7 ) 


where  d(i)  is  the  ith  data  point  of  the  above  par  ameters . 
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To  show  the  difference  between  mean-of-the-logarithms  and 
logarithm-of-the-mean,  (Section  II-G  and  Appendix  C),  the  value  LOG  MEAN 
was  computed  for  the  means  of  AMP,  RMS,  and  AMP/ RMS. 

Editing  of  the  noise  data  was  then  performed  by  rejecting  sam- 
ples for  which  at  least  one  of  the  values  AMP/ RMS,  LOG  AMP,  LOG  RMS,  or 
LOG  RATIO  fell  outside  an  interval  of  NSD  times  their  standard  deviation 
about  tie  mean.  The  number  NSD  can  be  any  number;  in  the  processing 
presented  here  NSD  = 2 was  selected.  This  editing  process  was  continued 
until  all  data  fell  within  the  desired  interval.  Final  editing  results  are  dis- 
played in  Appendix  E for  all  stations  involved. 

This  editing  process  may  not  be  stable.  If  the  data  contains 
many  high  extremes  it  may  force  out  the  low  values  and  bias  the  population 
high.  In  general,  noise  data  do  not  contain  anomalously  low  values,  but  the 
population  may  still  be  biased  low  if  the  rejected  high  values  were  established 
by  actual  seismic  noise.  If  the  distribution  of  values  is  approximately  normal, 
however,  the  above  criterion  is  fairly  reasonable.  Intentionally,  therefore, 
we  do  not  edit  the  maximum  amplitudes  or  RMS  values,  which  are  not  normally 
distributed,  but  rather  we  edit  their  logarithms  which  approximate  normal 
random  variables.  In  addition,  the  supposedly  near-constant  values  AMP/RMS 
were  edited  in  the  same  manner,  which  should  serve  especially  to  remove 
spikes  and  possibly  d.c.  levels. 

After  the  processing  was  completed  the  PDP-15  interactive 
display  at  the  SDAC  was  used  for  a quick  check  on  the  above  procedure  by 
means  of  a simple  visual  inspection  of  the  wide  band  noise  traces.  As  expect- 
ed, the  above  procedure  proved  by  no  means  perfect,  but  the  remaining  ano- 
malous data  and  the  rejection  of  some  presumably  good  data  were  judged  not 
to  significantly  distort  the  population  of  valid  noise  samples.  The  resulting 
data  displayed  in  Appendix  E,  at  least,  seem  plausible.  The  Z LP  waveforms 


appeared  somewhat  distorted.  The  V juPE  noise  study  by  Prahl  (1974)  dis- 
cusses the  general  editing  process  in  more  detail. 

C.  ESTIMATING  DETECTION  THRESHOLDS  FROM  NOISE 

The  detection  capability  of  a seismic  station  may  be  estimated 
by  deteri  lining  the  detection  thresholds  as  indicated  in  Section  II.  As  shown, 
we  need  not  necessarily  assume  a statistical  model  as  long  as  we  are  able  to 
determine  the  detection  probability  distribution  function,  F (x)  » otherwise. 
Since  we  want  to  compare  our  results  with  estimates  from  the  maximum  like- 
lihood direct  method,  which  assumes  a Gaussian  distribution,  we  also  assumed 

F (x)  to  be  normal.  As  pointed  out  in  the  analysis,  this  assumption  is  ap- 
M 

proximately  valid  if  the  range  of  relative  epicentral  distances,  /\/A,^  is  kept 
small.  The  available  maximum  likelihood  estimates,  however,  v'ere  deter- 
mined from  events  spread  over  relatively  large  regions,  and  for  most  stations 
the  log  \ distribution  is  not  normal,  not  even  where  the  event  populations 
were  divided  into  subsets  of  \<  50°  and  \ > 50°,  as  seen  from  Appendix  A. 
Since  both  methods  are  based  on  the  same  detection  mechanism  formulated  in 
Section  II,  they  will  make  approximately  the  same  error  in  making  the  assump- 
tion of  normality  for  the  probability  of  detection,  and  should  still  yield  similar 
results. 

As  was  shown  further  in  Section  II,  the  50%  detection  threshold 
is  determined  from  the  average  of  the  r.  v.  M.  This  means  that  we  must  use 
the  arithmetic  mean  of  the  logarithms,  or  equivalently,  the  logarithm  of  the 
geometric  mean  of  the  noise  amplitude,  AMP,  the  signal  period  T,  and  the 
epicentral  distance,  \ . We  furthermore  need  the  value  of  the  station  magni- 
tude bias  b,  and  an  estimate  of  the  s.  d.  for  the  parameters  log  AMP,  log  T, 
log  \ , and  for  the  station  M variation  or  bias  term  b . Furthermore,  the 
instrument  response  G(T),  and  the  constant  C = 1.  12  + detection  SNR  mar- 
gin are  required. 
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The  above  considerations  lead  to  the  following  procedure': 

• Determine  the  mean  and  s.  d.  of  the  logarithm  of  the  maximum 
peak-to-peak  seismometer  output  noise  amplitudes; 

MEAN  LOG  AMP,  LOG  AMP  S.  D. 

• Determine  from  previously  processed  events  in  the  region  of 
interest,  for  the  parameters  T,  \ , and  b ; 

the  logarithmic  dist ribution  of  T and  \ ; 

the  geometric  means  T , A ; 
n o o 

the  s.d.  of  the  logarithms;  O , <7  ; 

log  T log  A 

station  M bias  b and  s.  d.  (T 
s b 

• Input  to  the  detection  capability  estimator  the  values 

MEAN  LOG  AMP,  LOG  AMP  S.  D.  ; 

Tq  , G(To);  other  T , G(T)  values; 

G(20),  G ( i 0 ) ; 

, other  \'s  ; 

17  (BIAS)  and  O (BIAS  S.  D.  ) ; 
d(T)  = MS(20)  - MS(PER)  . 

The  detection  capability  estimation  algorithm  then  computes: 

• The  50%  detection  threshold. 

xcn  = = MEAN  LOG  AMP  - locT.G(T  ) 

50  M h o o' 

+ log  \ + d(T  ) - b + C ; (111-8) 

o o 

• The  standard  deviations: 

SIG  PER  EFFECTS  S.  D.  a]  = (/}-«-  1 ) a . _ ; (I1I-9) 

M log  i 

TOTALS.  D.  = a = [ (LOG  AMPS.D.)2 

M l 

+ (SIG  PER  EFFECTS  S.  D.  )Z  + (DIST  S.  D.  / 

+ (BIAS  S.  D.  )l  J 1/Z  , (111-10) 


111-9 


where 


1 


DIST  S.  D.  = (?,  * 

loe  A 


The  90%  detection  threshold: 


x = M +1.28(7  . 

90  M M 


(III- 11) 


The  value  MEAN  LOG  AMP  may  be  replaced  by  (MEAN  LOG 
RMS  + LOG  MEAN  RATIO)  if  the  use  of  RMS  values  rather  than  peak-to-peak 
amplitudes  is  preferred,  and  the  mean  ratio  is  known. 

Since  most  maximum  signal  amplitudes  were  found  to  have  a 
30-sec  period,  Figure  A-7,  and  since  detection  thresholds  are  given  as  20- 
sec  magnitudes,  we  computed  the  thresholds  for  20-sec  and  30-sec  signal 
periods,  rather  than  for  the  signal  period  geometric  mean.  The  results  then 
reflect  the  range  of  detection  capability  estimates  as  a function  of  signal  per- 
iod. Furthermore,  where  applicable,  we  computed  the  thresholds  for  the 
geometric  mean  of  the  A-values  for  the  total  region,  the  A <50  region, 
and  the  A > 50°  region.  Finally,  we  assumed  the  stations  to  be  unbiased  in 
magnitude  (STA  MS  BIAS  = 0.  0)  with  a s.  d.  of  0.  35  (BIAS  S.  D.  = 0.  35). 

D.  PRESENTATION  OF  THE  RESULTS 

The  results  of  the  editing  and  averaging  process  described  in 
Subsection  III-B  are  given  in  the  tables  of  Appendix  E.  Some  of  the  key  values 
from  these  tables  are  combined  in  an  overview  in  Table  III— 2 to  evaluate  the 
logarithm-of-the-mean  versus  the  mean-of-the-logarithm  operation  on  noise 
amplitudes,  and  the  consistency  of  AMP/  RMS  ratios.  The  threshold  esti- 
mates and  related  values  as  described  in  Subsection  III-C  are  listed  in  Tables 
HI- 3 through  III-8.  The  50%  and  90%  detection  threshold  estimates  from 
noise,  M50  and  M90,  are  given  for  the  case  that  noise  AMPLITUDE  values 
were  used  as  well  as  for  the  case  of  threshold  estimation  from  RMS  noise 
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TABLE  III  — 3 

DETECTION  CAPABILITY  ESTIMATION  FROM  NOISE:  STATION  CTA 


TABLE  III- 5 

DETECTION  CAPABILITY  ESTIMATION  FROM  NOISE:  STATION  KON 


I 


TABLE  III  - 7 

DETECTION  CAPABILITY  ESTIMATION  FROM  NOISE:  STATION  ALQ 


levels.  The  value  MT  represents  the  noise  magnitude  part  of  the  50% 
detection  threshold;  it  excludes  d(T)  and  the  bias  term  b: 

MT  = M = MEAN  LOG  AMP  - LOG  T G(T) 

N 

+ LOG  A + C.  (Ill- 12) 

o 

Furthermore,  the  tables  give  the  DIRECT  50%  and  90%  detection  threshold 
estimates  as  determined  by  the  maximum  likelihood  direct  method  from 
VLPE  data.  Table  III  - 9 gives  an  overview  of  50%  detection  threshold  differ- 
ences and  related  values. 

E.  DISCUSSION 

We  will  first  discuss  the  data  presented  in  Appendix  E ani 
in  Table  III -2.  The  log  amplitude  distributions  of  the  bandpassed  and  edited 
noise  data  are  given  in  Appendix  A.  They  are  approximately  normally  distri- 
buted and  their  standard  deviations  range  from  0.  1 to  0.  2 , Furthermore, 
we  notice  that  for  all  parameters 

| LOG  MEAN  - MEAN  LOG  | « 0.  2 1 5 ( <y//i  )2  (III- 13) 

as  predicted  in  Appendix  C.  For  most  stations  this  difference  is  less  than 
0.02  indicating  that  one  could  take  the  logarithm  of  the  arithmetic  rather 
than  the  geometric  mean  without  introducing  a significant  error.  This  also 
holds  for  the  AMP/RMS  ratio,  which,  however,  varies  more  than  expected 
at  the  station  CTA.  This  is  the  only  station  with  a definitely  non-log-normal 
noise  amplitude  distribution.  Thus,  we  may  estimate  detection  thresholds 
from  RMS  noise  levels  rather  than  from  noise  peak-to-peak  amplitudes  if 
preferred,  using  a fixed  AMP/PMS  ratio.  This  ratio  only  seems  to  be  con- 
sistent among  the  stations  KON,  KIP,  ALQ,  and  ZLP,  however,  and  the  value 
of  this  ratio  must  apparently  be  established  separately  for  each  station,  so 
*hat  this  method  offers  little  advantage,  especially  when  evaluating  a new  site. 
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Although  the  noise  amplitude  periods  (AMP  PER)  are  not  re- 
quired in  the  calculation  of  detection  thresholds  they  show  the  interesting 
fact  that  they  are  in  general  in  the  20  to  32  sec  range;  their  distributions  are 
given  in  Appendix  B. 

We  now  turn  to  the  results  presented  in  Tables  111-3  through 
111-9.  The  standard  deviations  of  the  combined  signal  period  effects  were 
already  given  in  Table  11-1;  together  with  the  log  amplitude  s.  d.  , the  log 
distance  s.  d.  , and  the  bias  s.  d.  , they  yield  a total  which  ranges  from 

0.39  to  0.9i,  in  accordance  with  Equation  (11-18). 

The  threshold  difference  due  to  using  either  20-sec  or  30-sec 
signal  periods  according  to  Equation  (A-6)  is  approximately: 

M90(2  0)  - M90(30)  « ()3-«-l)(log  20  - log  30) 

= -0.  18  (/3-a-l)  . (111-14) 

This  difference  is  greater  for  stations  with  signals  received  over  mixed 
continental-oceanic  propagation  paths  and  for  stations  with  a large  G(30)/G(20) 
instrument  response  ratio.  Over  the  six  stations  used  in  the  evaluation,  the 
difference  ranges  from  0.  18  (KON)  to  0.  98  (CT A)  as  is  also  reflected  by 
Table  A-l.  Since  the  station  magnitude  bias  is  assumed  zero,  the  detection 
threshold  for  a 20-sec  signal  period  is  given  by  the  noise  threshold  magnitude 
MT  ; for  30-sec  signal  period  detections  we  must  add  the  period  dependent 
magnitude  difference  terms  given  by  Equation  (111-14). 

According  to  the  overview  in  Table  111  - 9 , the  90%  detection 
thresholds  show  excellent  agreement  between  the  estimates  obtained  from 
noise  samples  and  the  maximum  likelihood  direct  estimates.  For  three  of 
the  six  stations  (CHG,  KON,  KIP)  the  direct  estimate  lies  about  midway  be- 
tween the  20-  and  30-sec  estimates  from  noise.  According  to  the  confidence 
level  in  the  maximum  likelihood  detectability  curves  (Lambert  et  al.  , lp7  3) 
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the  direct  estimates  for  these  stations  seem  reliable.  This  may  indicate  that 
the  geometric  mean  of  the  periods  of  maximum  signal  amplitude  for  those 
stations  is  approximately  25  seconds. 

The  CT A estimate  from  noise  seems  somewhat  high,  although 
the  direct  estimate  still  lies  between  the  20-  and  30-sec  signal  period  noise 
thresholds.  This  may  be  related  to  the  strong  total  effect  of  the  signal  per- 
iod due  to  the  combination  of  a steep  instrument  response  and  the  fact  that 
the  station  receives  its  signals  over  mixed  continental-oceanic  paths.  Also, 
the  nc  e amplitude  distribution  is  not  log-normal,  and  the  log  amplitude  s.  d. 
is  higher  than  those  of  CHG,  KON,  KIP,  ALQ,  and  7,LP.  The  direct  esti- 
mates for  this  station  may  be  less  reliable  due  to  a relatively  low  number  of 
higher  magnitude  events;  the  confidence  limits  seem  somewhat  wide. 

The  stations  ALQ  and  7 LP  show  a 0.  1 to  0.  3 lower  detection 
threshold  estimated  from  noise  as  compared  with  the  direct  estimate.  ALQ 
noise  traces  looked  good  on  the  interactive  display;  its  log  amplitude  distri- 
bution is  close  to  normal,  and  its  log  amplitude  s.  d.  and  the  combined  signal 
period  effects  s.  d.  are  small.  The  ALQ  direct  estimate  may  be  less  reliable 
than  that  of  the  other  stations,  since  the  population  of  processed  seismic 
events  seems  to  contain  relatively  few  high  magnitude  events;  the  confidence 
limits  are  wider  than  those  for  CTA.  The  same  situation  exists  to  a greater 
extent  for  the  station  ZLP,  for  which  the  overall  population  of  processed 
events  is  low.  Also,  the  display  of  noise  traces  consistently  showed  a rather 
unusual  waveform  for  this  station,  which  may  affect  the  detection  capability 
estimates  from  noise  as  well  as  the  direct  estimates.  Nevertheless,  the  log 
amplitude  distribution  is  close  to  normal  and  since  the  waveform  shape  may 
not  disturb  the  detection  mechanism  very  much,  the  estimates  from  noise 
may  well  be  reliable. 

The  90%  detection  thresholds  are  difficult  to  compare  since, 
if  the  50%  thresholds  agree,  the  90%  threshold  differences  depend  solely  on 
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the  standard  deviation,  a , used  in  x = x&Q  + 1.28  a . As  discussed  earlier, 
the  standard  deviation  used  in  maximum  likelihood  direct  estimation  obtained 
via  regression  on  the  bodywave  to  surface-wave  magnitude  relationship  seems 
low  with  respect  to  the  s.d.  estimated  from  the  statistics  of  each  term  in- 
volved in  estimating  detection  thresholds  from  noise.  This  is  reflected  in  the 
results  given  in  Tables  III- 3 through  III-8  where  the  TOTAL  S.  D.  in  estima- 
tion from  noise  should  be  compared  with  the  value  a = 0.  4 used  in  the  direct 
method. 

From  the  above  considerations  we  may  conclude  that  estima- 
ting the  detection  capability  from  ambient  noise  is  a valid,  fast,  and  accurate 
method.  A population  of  approximately  30  valid  noise  samples  already  gives 
reliable  estimates.  This  means  that  for  a well  performing  station,  taking 
noise  samples  probably  once  or  twice  a day,  and  accounting  for  50%  rejection 
(e.  g.  , KIP),  a reliable  detection  capability  estimate  can  be  established  within 
one  to  two  months. 

The  method  also  required,  however,  statistics  on  the  periods 
of  maximum  event  signal  amplitudes,  the  epicentral  distances,  and  the  station 
magnitude  bias.  For  existing  stations  these  statistics  may  be  obtained  from 
previously  processed  events.  For  a new  or  proposed  site  the  epicentral  dis- 
tances can  be  calculated  from  the  seismicity  in  the  region  of  interest,  where- 
as the  signal  period  and  magnitude  bias  statistics  can  only  be  estimated,  as- 
suming that  these  statistics  will  be  similar  to  those  of  the  existing  stations. 
Among  the  VLPE  stations  the  signal  period  statistics  looked  quite  similar; 
there  is  yet  insufficient  information  about  the  station  magnitude  bias  statistics. 
In  the  evaluation  presented  here  the  signal  period  statistics  of  the  station  KIP 
were  considered  typical,  and  were  applied  to  ail  six  stations.  Furthermore, 
a zero  magnitude  bias  with  a s.  d.  of  0.  35  was  assumed  at  each  station.  Evi- 
dently, that  procedure  worked  well  for  ail  stations  used  in  the  evaluation  and, 
unless  the  physics  (path,  local  structure,  instrumentation)  of  a new  station 
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are  definitely  different,  it  is  expected  that  estimating  signal  period  and  mag- 
nitude bias  statistics  from  the  statistics  of  existing  stations  will  not  introduce 
a significant  error  in  estimating  the  new  station's  detection  capability  from 
noise. 

Finally,  it  should  be  noted  that  the  accuracy  of  the  detection 
thresholds  estimated  from  noise  was  established  with  respect  to  the  maxi- 
mum likelihood  estimates.  As  pointed  out  previously,  both,  the  maximum 
likelihood  method,  and  the  estimation-from-noise  procedure  presented  here 
assume  a Gaussian  model  for  the  probability  of  detection.  A deviation  from 
this  model  by  any  of  the  random  variables  involved  would  cause  approximate- 
ly the  same  error  in  either  method.  Thus,  although  these  two  methods  show 
ex  ellent  agreement,  they  may  both  be  in  error  with  respect  to  the  true 
probability  of  detection. 


SECTION  IV 


SUMMARY 


If  detection  of  a seismic  event  at  a given  station  is  based  on 
the  surface  wave  signal-to-noise  ratio  (SNR)  the  probability  of  detecting  an 
event  given  its  20-sec  reference  magnitude  is  given  by  the  probability  dis- 
tribution function  of  a r.  v.  M : 


P(det  M (20)  = x)  = P(M<x)  = F (x)  , 
S IVl 


(IV-1) 


where 


M (20)  is  the  20-sec  reference  magnitude  for  the  event  con- 
s 


cer ned; 

is  the  numerical  value  of  the  reference  magnitude; 


F (x)  is  the  probability  distribution  function  of  the  r.  v,  M . 
M 


The  r.  v.  M is  described  by 


M = Mn  + d(T)  - b , 


(IV  -2) 


where 


M 


N 


is  the  so-called  "noise  magnitude"  ; 
d(T)  is  the  20-sec  minus  T-sec  magnitude  difference; 
T 


is  the  period  of  the  maximum  event  signal  amplitude  at 
the  seismometer  output  ; 


is  the  station  magnitude  variation  about  the  ».  eference 
magnitude. 


IV-1 


i 


3 


I 


r 


The  "noise  magnitude"  is  given  by 


T 


where 


= l°g  ANout  " lo8  T.G(T)  + log  A + C , 


(IV-3) 


Am  is  the  maximum  peak-to-peak  noise  amplitude  at  the 
lNout 

seismometer  output; 

T is  the  signal  period  defined  above; 

G(T)  is  the  instrument  response  for  the  period  T , scaled 

so  that  a 40- sec  output  amplitude  represents  true 
ground  motion  in  millimicrons,  G(40)  = 1 ; 

A is  the  epicentral  distance  for  the  event  and  station  con- 

cerned, 

C = 1.12  4 a SNR  detection  margin  (e.g.  , 0.15  for  a 

3dB  rrargin.) 


The  term  "noise  magnitude"  is  not  quite  appropriate  since  it  consists  of  a 

no' se  parameter  (Am  ) on  the  one  hand,  and  event  signal  parameters  (*,A) 
r tNout 

on  the  other. 


The  above  model  of  the  detection  mechanism  does  not  presume 
any  special  statistical  characteristics,  and  therefore  may  serve  as  a model 
for  any  method  of  estimating  a seismic  station's  detection  capability.  Based 
on  this  model  two  methods  of  estimating  detection  capabilities,  and  two  ap- 
proaches to  each  method  were  discussed: 

• Estimation  from  actual  detection  percentages 

- not  assuming  statistical  characteristics  (the  so-called 
direct  method) 

- assuming  a Gaussian  probability  of  detection  (the  maxi- 
mum likelihood  approach  to  the  direct  method). 
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Estimation  from  noise 


not  assuming  specific  statistical  characteristics; 
assuming  a Gaussian  probability  of  detection. 

Available  data  indicates  that  the  probability  of  detecting  an 
event  given  its  magnitude  can  be  considered  to  be  approximately  Gaussian,  if 
the  relative  distance  range,  A/Aq  , is  kept  small  (e.g.  , A = 50°, 

35° < A < 80°  ; A = 100°  , 70°<  A < 160°). 

o ' 

In  the  general  case  the  50%,  90%,  etc.,  detection  thresholds 
(X50*  xgo*  etc*  ) are  f°und  from 


fm<V 


= 0.5 


and 


FM(x90> 


= 0.  9 . 


(IV-4) 


(IV-5) 


The  Gaussian  probability  of  detection  is  described  by 


Fm(x>  =0.5+  erf 


x 


(IV-6) 


where  and  are  the  mean  and  standard  deviation,  respectively,  of 

the  r.  v.  M.  In  this  case  the  detection  thresholds  are  determined  by  M and 


and 


50 


= V 


M 


90 


|i..  + 1.28 
M 


M 


(IV-7) 


(IV-8) 


In  estimating  detection  capabilities  from  noise  the  statistics  of 
each  term  constituting  the  r.  v.  M,  i.  e.  , the  statistics  of  noise  amplitude, 
event  signal  period,  event  epicentral  distance,  and  station  magnitude  bias, 
are  combined  to  yield  the  detection  probability  distribution  function  F (x) 
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In  the  general  case  the  entire  F fx)  curve  must  be  established;  in  the  Gaus- 
° M 

sian  case  we  only  need  to  determine  the  mean  and  s.  d.  of  M . 

Instrument  output  RMS  noise  values  may  be  used  instead  of 
maximum  noise  amplitudes  in  this  method  if  the  average  maximum-amplitude- 
over-RMS  ratio  is  known  for  the  station  concerned.  This  ratio  was  approxi- 
mately equal  among  four  out  of  six  VLPE  stations. 

Based  on  the  analysis  of  the  two  methods  and  their  approaches 
it  was  concluded  that  for  the  non-Gaussian  as  well  as  for  the  Gaussian  case 
there  should  be  little  difference  between  the  estimates  obtained  from  actual 
detection  percentages  and  those  obtained  from  noise  samples.  This  conclu- 
sion was  supported  by  an  empirical  evaluation  performed  on  six  stations  of 
the  V LPE  network,  by  comparing  detection  capabilitie  s estimated  from  noise, 
assuming  a Gaussian  detection  model,  with  estimates  obtained  by  the  maxi- 
mum likelihood  method.  The  50%  detection  threshold  estimates  of  the  two 
methods  agree  in  general  within  0.2  magnitude  units.  The  90%  detection 
thresholds  differ  due  to  a difference  in  variance  between  the  two  methods, 

since  the  maximum  likelihood  variance  was  obtained  through  m,  to  M con- 

b s 

version.  The  minimum  statistical  population  for  estimating  detection  capa- 
bilities from  noise  appears  to  be  in  the  order  of  30  valid  noise  samples,  and 
a station's  detection  capability  may  be  established  within  one  or  two  months 
for  a well-performing  instrument. 

This  technique,  however,  requires  a prior  knowledge  of  the 
statistics  of  the  periods  of  maximum  signal  amplitudes,  the  epicentral  dis- 
tances, and  the  station  magnitude  bias.  For  an  existing  station  these  statis- 
tics may  be  obtained  from  previously  processed  events.  For  a new  or  pro- 
posed site  the  distance  statistics  can  be  established  from  the  seismicity  in 
the  region  of  interest.  The  other  statistics  may  be  estimated  from  the  event 
processing  performed  for  other  stations  of  a similar  physical  configuration. 
Furthermore,  the  estimates  obtained  from  noise  may  require  seasonal  ad- 
justments. 
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Finally,  it  is  pointed  out  that  deviation  of  the  statistics  con- 
cerned from  the  Gaussian  detection  model  will  introduce  similar  errors  in 
the  maximum  likelihood  estimates  and  in  the  estimates  obtained  from  noise. 
Agreement  between  these  two  methods,  therefore,  does  not  necessarily 

guarantee  that  these  estimates  are  correct  with  respect  to  the  actual  detec- 
tion capabilities. 

Related  topics  that  need  further  investigation  are  the  distribu- 
tions and  mutual  correlations  of  sigr.*!  period,  epicentral  distance,  and  sta- 
tion surface  wave  magnitude  bias;proCedures  for  the  selection  of  valid  noise 
samples;  extension  of  this  study  to  bodywave  magnitudes  and  the  surface- 

wave  horizontal  components;  and  comparison  of  non-Gaussian  detection 

capability  estimation  from  noise  with  the  non-Gaussian  detection  estimates 
of  the  diret  t method. 
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APPENDIX  A 


STATISTICS  OF  THE  RANDOM  VARIABLE  M 


This  appendix  investigates  the  statistics  of  the  terms  and  their 
combinations  contributing  to  the  r.  v.  M given  by  Equation  (11-16): 


where 


M = log  A^  ^ - log  T - log  G(T)  + log  A + C + d(T)  - b (A-l) 
is  the  maximum  peak-to-peak  noise  amplitude  at  the 

1Nout 

seismometer  output; 

T is  the  period  cf  the  maximum  event  signal  amplitude; 

G(T)  is  the  instrument  response  for  the  period  T ; 

A is  the  event  epicentral  distance; 

C combines  the  distance  correction  term  (1.  12)  and  a de- 

tection margin  term; 

d(T)  is  the  20-sec  minus  T-sec  magnitude  difference; 

b is  the  station  magnitude  variation  with  respect  to  the 

reference  magnitude. 


The  following  subsections  discuss  the  statistics  of  the  terms:  log  ^Nout  » 
log  T,  log  G(T),  d(T),  the  combination  d(T)  - log  T - log  G(T),  and  the  terms 
log  A and  b , respectively.  The  final  subsection  combines  this  information 
to  an  indication  about  the  statistics  of  the  r.  v.  M . Indications  about  the  sta- 
tistical distribution  of  most  of  these  terms  are  based  on  visual  inspection  of 
histograms  given  for  the  parameters  log  ^•Nout-’  T,  anc*  l°g  A . 


A-l 


1. 


Statistics  of  the  Maximum  Peak-to-Peak  Noise  Output  Ampli- 
tude Ai\t 

Nout 

To  obtain  an  indication  of  the  statistical  distribution  of  log  am- 
plitude values  for  he  maximum  noise  amplitude  at  the  seismometer  output, 
the  histograms  of  the  log  amplitudes  actually  used  in  the  empirical  evaluation 
Part  of  this  study  are  given  in  Figures  A- I through  A-6.  These  maximum 
amplitudes  were  obtained  from  a population  of  noise  samples  in  the  passband 
17-44  sec.  Die  population  was  edited  lo  remove  extreme  values  as  described 
in  Section  III.  The  maximum  amplitude  distributions  for  the  VLPE  stations 
CHG,  KON(  KIP,  ALO,  and  ZTP  may  be  considered  to  be  approximately 
log-normal.  The  reason  for  the  somewhat  different  amplitude  behavior  at 
the  station  CTA  is  not  known  at  this  point,  but  may  be  partially  caused  by 
insufficient  editing. 

Another  possible  indication  of  the  log-normality  of  maximum 
noise  amplitudes  stems  from  the  log-normality  of  RMS  noise  levels  at  the  in- 
strument output,  observed  by  Alsup  and  Becker  (1973).  In  Appendix  D it  is 
shown  that,  if  certain  assumptions  arc  proven  valid,  the  log-normality  of 
RMS  input  noise  levels  implies  the  log-normality  of  output  maximum  noise 
amplitudes,  using  the  fact  that  the  logarithm  of  most  long  period  instrument 
responses  is  either  nearly  linear  with  the  period  logarithm  (VLPE  stations), 
or  almost  constant  (ALP A,  NORSAR)  in  the  10-  to  30-sec  period  range. 

Thus,  the  available  data  indicate  that  in  general  the  maximum 
noise  amplitude  at  the  seismometer  output  can  be  considered  to  have  an  ap- 
proximately log-normal  probability  distribution. 

2.  The  Signal  Period  T 

'1  o determine  the  influence  of  the  period,  T , of  the  maximum 
signal  amplitude,  histograms  were  compiled  from  VLPE  detection  data  over 
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the  year  1972.  A typical  histogram,  the  one  for  station  KIP,  given  in  Figure 
A-7,  shows  that  most  maximum  signal  amplitudes  have  a period  of  30  seconds, 
and  smaller  but  relatively  significant  numbers  of  maximum  amplitudes  are 
found  at  20,  22,  and  28  seconds.  This  results  in  the  log  T histogram  of  Fig- 
ure A-8,  which  indicates  that  the  signal  period  distribution  may  be  consider- 
ed to  be  approximately  log-normal.  The  standard  deviation  of  log  T appears 
to  be  in  the  order  of  0.  1 for  most  V LPE  stations. 

3.  The  Instrument  Response  G(T) 

Fjr  most  V LPE  stations  the  logarithm  of  the  instrument  re- 
sponse is  approximately  linear  with  the  logarithm  of  the  period  in  the  20-  to 
30-sec  period  range  (Figure  A-9): 

log  G(T)  =:  or  log  + log  G(20),  20  < T < 30  , (A-2) 

where 

„ log  G( 30)  - log  G(20)  log  G(30)/G(20) 

log  30  - log  20  0.  18  ' ' 

Since  a linear  transformation  of  a normal  r.  v.  is  again  a normal  r.  v.  (e.  g.  , 

Papoulis,  1969),  and  log  T is  approximately  normal,  the  r.  v.  log  G(T)  is 
also  approximately  normal. 

The  long  period  arrays  of  ALPA  and  NORSAR  have  an  almost 
flat  response  between  20  and  30  seconds,  so  that  tnere  log  G(T)  is  approxi- 
mately constant  and  does  not  affect  the  statistics  of  the  r.  v.  M. 

4.  Statistics  of  the  Magnitude  Difference  Due  to  Period  Effects, 
d(T) 

The  difference  due  to  period  and  path  effects  combined  is 
plotted  in  Figure  A- 10,  after  data  provided  by  Marshall  and  Basham  (1972). 

The  period  is  plotted  on  a logarithmic  scale.  The  curves  appear  to  be  approxi- 
mated' linear  in  segments  between  20  and  30  seconds  and  30  and  49  seconds. 
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Since  according  to  Figure  A-7,  most  signal  periods  fall  in  the  20-  to  30-sec 
range,  and  since  the  signal  period  has  an  approximately  log-normal  distribu- 
tion, the  magnitude  difference  will  also  be  approximately  normal.  The  near- 
linearity with  log  T of  the  magnitude  difference  d(T)  due  to  period  effects 
only,  may  be  expressed  by 

d(T)  ~ ft  log  -J5  * 20  < T £ 30.  (A-4) 

For  Eurasian  continental  paths  ft = 1.70,  for  mixed  oceanic -continental  paths 
ft  = 0.  68  according  to  the  Marshall  and  Basham  data. 

The  path  effects  shown  in  these  curves  presumably  are  absorb- 
ed in  the  station  magnitude  bias  b . 

5.  Statistics  of  the  Combined  Signal  Period  Effects  on  M 

It  now  is  of  interest  to  study  the  total  effects  on  the  r.  v. 

M = + d(T)  - b due  to  the  statistics  of  the  signal  period  T.  In  particular, 

we  are  interested  in  the  difference  between  the  T-sec  and  20-sec  M values 
determined  by  the  period-dependent  terms  log  T and  log  G(T)  in 
(Equation  11-10),  and  by  the  term  d(T)  in  Equation  (11-15): 

M(T)  - M(2 0)  ~ d(T)  - log  ~ - log  , 20<T<30.  (A-5) 

thus,  the  empirically  determined  magnitude  difference  due  to  period  effects, 
d(T),  is  counteracted  by  the  terms  log  T and  log  G(T). 

With  Equations  (A-3)  and  (A-4)  the  above  relationship  becomes: 

M(T)  - M(20)  ~ (/3  - or  - 1 ) log  — , 20<T<30.  (A-6) 

Since  the  combined  period  effects  are  approximately  linear  with 
log  T the  standard  deviation  of  the  combined  period  effects  can  be  expressed 
in  terms  of  the  standard  deviation  of  log  T : 


A-  1 0 


I . 


L 1 


°M~  l'J",-1lCTlos  T • (A-7) 

whe  re 

T 

CTM  is  the  standard  deviation  of  the  combined  signal  period 
dependent  terms  in  M ; 


* T *S  t*ie  stauc*arc*  deviation  of  log  T . 

T 

Thus  the  value  of  <7  depends  on  the  signal  period  distribution,  on  the  slope 
of  the  instrument  response  between  20  and  30  seconds,  and  on  the  rate  of  mag- 
nitude change  with  the  signal  period  logarithm.  For  a flat  instrument  response 
(ALPA,  NORSAR)  a equals  zero  in  the  above  equations. 


Table  A- 1 lists  the  20-  and  30-sec  instrument  responses  for  the 

period  November  1972  through  January  1973,  and  the  values  n , (3  , fi-a-l  , 

T 

anC*  °M  ^°r  ^ LPE  stations  used  in  the  empirical  evaluation  ol  this  study 

using  a (X^  T velue  of  0.  1 (Figure  A- 8 ).  The  stations  CIIG,  KON,  ALQ, 
and  ZLP  thereby  are  assumed  to  receive  most  Eurasian  event  signals  over 
continental  propagation  paths;  the  paths  to  the  stations  CT A and  KIP  are  in 
general  mixed  continental-oceanic.  We  notice  that  the  log  T G(T)  term  al- 
ways overcompensates  the  term  d(T)  ; the  extent  of  overcompensation  de- 
pends on  the  factors  mentioned  above. 


6.  The  Epicentral  Distance  \ 

We  now  turn  to  the  next  variable  in  the  noise  magnitude,  the  epi- 
central distance  \ . For  a relatively  large  area  we  would  not  expect  the  para- 
meter A to  have  a log-normal  distribution.  For  example,  for  the  population 
of  events  (detected  and  undetected)  processed  by  Lambert  et  al.  (1973)  the 
log  \ distribution  of  the  VLPE  stations  is  given  in  Figures  A- 11  through 
A-22.  For  instance,  for  KIP  (Figure  A-18),  the  distances  range  from  14°  to 
146°,  and  their  logarithms  accordingly  from  1.4  to  2.2.  In  processing  the 
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events  Lambert  ei  al.  (197  3)  split  the  \ -populations  into  groups  of  events 
with  epicentral  distances  either  greater  or  less  than  50°.  For  KIP  there 
distributions  are  given  in  Figure  A-19  and  Figure  A-20,  respectively.  1 he 
bin  size  causes  an  overlap  betv/een  the  two  groups  at  log  SO  = 1.7  . 

It  is  noticed  that  decreasing  the  \ - range  in  the  above  man- 
ner leads  only  in  a few  cases  to  approximation  of  a log-normal  distribution. 
Thus,  if  a log-normal  A -distribution  is  desired  for  a certain  detectability 
estimation  model,  a relatively  narrow  distance  range  must  be  considered. 

In  the  case  of  direct  detect^oility  estimation  this  -nay  cause  the  population  of 
events  to  drop  below  the  statistically  required  minimum.  When  estimating 
the  detectability  from  noise,  however,  we  may  compute  detectability  thresh- 
olds for  one  particular  distance  value  A or  for  a relatively  narrow  \- 

o 

range  centered  about  some  A , and  shift  the  thresholds  over  an  amount 

log  \ /A  for  other  \ -values  of  interest.  The  variance  of  M due  to  distance 
o 

.variation,  of  course  increases  with  the  range  of  distances,  thereby  changing 
the  shape  of  the  detectability  curve  P(det  M^(20)  = x)  = F^(x). 

7.  The  Station  Magnitude  Variation  b 

The  station  magnitude  variation,  b,  about  the  reference  magni- 
tude is  described  by  Equation  (Il-h),  The  mean  variation  is  called  thr  station 
magnitude  bias.  Bodywave  magnitudes  were  found  to  be  approximately  nor- 
mally distributed  (Freedman,  1967).  For  detectability  estimation  modeling 
purposes  we  will  assume  the  variation  of  the  surface-wave  magnituces  also  to 
be  normal;  however,  this  needs  further  confirmation. 

8.  Total  Statistics  of  the  Random  Variable  M = M ^ + d('T')  - b 

With  the  statistics  known  on  each  of  the  terms  involved  in  the 
r.  v.  M we  are  now  able  to  determine  the  total  statistics  ot  this  variable.  We 
found  that  each  term,  except  the  epicentral  distance,  was  approximately  nor- 
mally distributed,  while  the  epicentral  distance  would  either  be  kept  fixed,  or 
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allowed  to  be  taken  only  over  a narrow  range  to  approximate  a log-normal 
distribution.  The  linear  combination  M of  these  approximately  normally 
distributed  random  variables  then  again  establishes  approximately  a normal 
r.  v.  . 

The  next  point  is  to  check  if  the  parameters  Aj^  T , \ , 
and  b are  correlated.  The  maximum  noise  amplitude  at  the  seismometer 
output,  Ajyj  is  not  correlated  with  any  of  the  parameters  T,  A , or  b.  The 
signal  period  T may  be  slightly  correlated  with  A because  of  path  effects. 
For  long-period  signals  and  small  event  areas,  however,  this  correlation  is 
expected  to  be  low.  There  may  also  be  some  correlation  between  T and  b, 
since  the  station  magnitude  variation  is  caused  by  radiation  pattern  and  pro- 
pagation path  effects.  Also  this  correlation  is  expected  to  be  low,  however. 
Finally,  for  the  same  reason,  there  may  be  some  correlation  between  A 
and  b . A thorough  study  of  these  correlations  would  be  desirable  but  we 
suggest  that  this  be  done  as  a separate  topic. 

If  all  variables  are  considered  to  be  mutually  uncorrelated, 
the  standard  deviation  of  the  r.  v.  M can  be  computed  from 


= + (/3-a-l)2  + a2  + a2  , 

M log  A log  T log  A b 


(A- 8) 


where 


or  is  the  standard  deviation  of  the  r.v.  M ; 

M 

C<  . is  the  standard  deviation  of  log  Am  : 

log  A 6 1Nout 

cr,  _ is  the  standard  deviation  of  log  T ; 

log  T K 

cr  . is  the  standard  deviation  of  log  A ; 

log  A 

cr  is  the  station  magnitude  bias  standard  deviation. 


The  standard  deviation  of  tho  logarithm  of  the  maximum  noise 
amplitude  at  the  seismometer  output,  cr^  must  be  determined  empirically 
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from  noise  sample  populations.  For  well-performing  VLPE  stations,  this 
value  turns  out  to  be  0.  1 tp  0.  2 (Section  III  and  Appendix  E).  The  stand- 
ard deviation  of  the  signal  period  logarithm,  CT^  y , was  found  to  be  approxi- 
mately 0.  1.  The  standard  deviation  of  the  combined  period  effects  then  fol- 
lows from  Equation  (A- 7);  their  values  were  listed  in  Table  A-l,  and  range 

from  0.  10  to  0.  33  depending  on  the  kind  of  path,  on  the  slope  of  the  instrument 
response,  and  on  the  slope  of  the  variable  d(T). 

The  standard  deviation  of  the  event  epicentral  distance  loga- 
rithm depends  on  the  relative  distance  range,  i.  e.  , it  depends  on  the  size  of 
the  area  of  interest  and  on  the  distance  between  that  area  and  the  station  con- 
cerned. For  most  VLPE  stations  and  Eurasian  events  a,  . ranges  from 
0.  to  0.2;  the  log  \ distribution  is  non-normal  for  such  a large  relative 
distance  range.  For  a fixed  A -value,  naturally,  cr^  0;  U is  small  for 
a narrow  relative  distance  range.  In  the  empirical  evaluation  of  this  study 
the  \ ranges  were  chosen  the  same  as  in  the  maximum  likelihood  direct 
detectability  estimation  reported  bv  Lambert  et  al.  (1973).  The  log  A dis- 
tributions for  these  cases  were  given  in  Figures  A-ll  through  A-22.  The 
station  magnitude  standard  deviation,  finally,  according  to  this  same  reference, 
v as  calculated  at  approximately  0.  35. 


A , and  b 


Assuming  no  correlation  among  the  parameters  , AM  T 

^out  9 

the  standard  deviation  of  the  r.  v.  M then  is  in  the  order  of 


o r 


iT 


M 


~ £(1.  0 to  0.  2)2  + (0.10  to  0.33) 


2 2 
+ (0.  1 to  0.  2) 


+ 0.  35 


2 J 1/2 


( A-9) 


0.  39  < O < 0.55. 
M ~ 


(A-l  0) 
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The  above  is  only  a preliminary  determination  of  the  statistics 
involved.  More  detailed  studies  should  be  performed,  both  analytically  and 
empirically,  to  further  test  the  assumptions,  observations,  and  calculations. 
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APPENDIX  B 

VLPE  LOG  NOISE  PERIOD  DISTRIBUTION 


The  following  pages  contain  the  distributions  of  the  logarithms 
of  the  maximum  noise  amplitude  periods. 
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APPENDIX  C 

SOME  CHARACTERISTICS  OF  THE  LOGARITHM 


Below  we  will  develop  some  characteristics  of  the  logarithm 
which  are  of  interest  in  the  computation  of  magnitudes  of  seismic  events  and 
their  averages. 

1.  Variation  of  Log  x Due  To  Variation  Of  x 

Consider  the  variable: 


then 


X = X + A X = 
o 


log  X = log  x^  + 


X ( 1 + 

o 

\ X 

X 

o 

log  ( 1 + 

\ X 

) . 


(C-l) 


(C-2) 


Thus,  the  variation  in  log  x , due  to  the  variation  in  x , is  determined  by 
the  relative  change,  \ k/x^  , and  equals 


\ x 


A(log  x)  = log  ( 1 + — — ) 


(C-3) 


For 


\x 


< 1 the  right-hand  side  may  be  developed  into  a Taylor  series,  so 


that  Equation  (C-2)  becomes 


log  x = log  x^  + 0.  4 


r \* 


c-i 


l 

2 ! 


1 / \x\  ^ 

~ "'■■-I  ' 


(C-4) 


For  small  values  of  the  variation  of  log  x with  x is  approximately 

o 

linear : 


log  x =:  log  x^  + 0.4  3 — 


\ x 


(C-5) 


For  a value  ---  <r  0.  5 , for  instance,  the  error  is  less  than  0.  OS.  The 
linear  approximation  and  the  exact  function  of  the  term  log  ( 1 + ) are 

plotted  in  Figure  C-  l. 

2.  Mean-of-the-Logarithm  Versus  Logarithm-of-the-Mean 

Assume  a r.  v.  x to  vary  symmetrically  about  its  mean.ji, 
with  a standard  deviation  cr  : 

\ * t (C-6) 


x = M Ax  = M(1  + 


) • 


Applying  Equation  (C-4)  results  in 

\ x 1 / _A_ 

log  x = log  +0.43  —jj  2 ! V M 


- (v)  • | 


(C  -7 ) 


To  obtain  the  mean  of  the  logarithm  we  take  the  expected  value  on  both  sides: 

2 


^ = log  ^ + 0.  4 3 E 

log  x 


Ax  _L 

— ■ 2 ! \ 4 ) 


: (V)  ••••;. 


(C -6) 


2 2 


Since  for  a symmetric 
this  becomes 


ic  distribution  of  \x  about  ^ E(\x)-  0 and  E(\x)  fT 


It  = log  It  - 0.215 

log  x 


(*)' 


+ higher  order  terms.  ((  -0) 


C-2 


FIGURE  C-l 

VARIATION  IN  LOG  x , DUE  TO  VARIATION  IN  x 


Thus,  for  small  cr//i  values  the  arithmetic  mean  of  the  logarithm  approxi- 
mately equals  the  logarithm  of  the  arithmetic  mean;  the  error  term  is  ap- 


p  roximately 


(C-10) 


For  instance,  for  a ratio  of  (7//i  = 0.5  the  error  term  is  approximately  0.05. 


The  validity  of  the  above  approximation,  and  in  particular  the 
accuracy  of  the  error  term,  depend  on  the  rate  of  convergence  of  the  higher 
order  terms,  which  is  a matter  of  third,  fourth,  etc.  statistical  moments. 
According  to  the  results  presented  in  Table  III  - 2.  the  approximation  of  the 
error  term  evidently  is  very  good. 
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APPENDIX  D 
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DISTRIBUTION  OF  INPUT  VERSUS  OUTPUT  NOISE  AMPLITUDES 


Alsup  and  Becker  (1973)  found  that  the  RMS  true  ground  motion 
noise  levels  were  approximately  log-normally  distributed.  Lambert  et  aJ. 

(197  3),  observed  that  at  the  seismometer  output  the  maximum-noise-amplitude- 
over-RMS  ratio  was  relatively  constant  for  V LPE  data.  We  will  show  that  if 
this  ratio  is  also  constant  at  the  instrument  input,  then  the  output  RMS  noise 
levels  and  maximum  noise  amplitudes  are  also  approximately  log-normally 
distributed. 


by: 


The  input  and  output  amplitudes  and  RMS  values  are  related 


log  Am  = log  RMS.  + log  (Am.  /RMS.  ) 
1Nout  in  ‘in  in 


+ logG(TN)  , (D-l) 

whe  re 

Ant  is  the  maximum  peak-to-peak  noise  output  amplitude; 

nOllt 

An.(  is  the  maximum  peak-to-peak  noise  input  amplitude; 

RMS,  is  the  RMS  input  noise; 
in 

G(T  ) is  the  instrument  response  for  the  period  T.  of  ‘he 
maximum  noise  amplitude. 

According  tc  the  tables  in  Appendix  E,  and  the  histograms  in 
Appendix  B,  the  noise  periods  T^  fall  mainly  within  the.  20-  to  30-see  range 
and  can  be  considered  to  be  approximately  log-normally  distributed.  Since 


D-l 


P 


4 


for  the  20-  to  30-sec  range  log  G(TN)  is  approximately  linear  with  log  TN> 
also  G(Tn)  will  be  approximately  log-normally  distributed.  If  the  input 
amplitude-over-RMS-noise  ratio  is  constant,  the  right-hand  side  of  Equation 
(D-l)  is  a linear  combination  of  two  normally  distributed  random  variables 
plus  a constant,  so  that  log  Atf^  must  be  approximately  normally  distrib- 
uted. Since,  furthermore, 

log  ANQut  = 108  RMSout  + 108  <ANout/RMSout)  * (D‘2) 

and  the  Ant  /RMS  ratio  was  found  to  be  relatively  constant,  the  out- 
Nout  out 

put  RMS  noise  values  also  arc  log-normally  distributed. 
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TABLE  E-4 

INTERMEDIATE  VALUES  USED  IN  DETECTION  CAPABILITY  ESTIMATION 

FROM  NOISE:  STATION  KIP 
(PAGE  1 OF  2) 
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